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Abstract 

In a /c-dimensional system of weakly interacting Bose atoms trapped by a 
spherically symmetric and harmonic external potential, an exact expression is 
obtained for the rotating ground states at a fixed angular momentum. The re- 
sult is valid for arbitrary interactions obeying minimal physical requirements. 
Depending on the sign of a modified scattering length, it reduces to either 
a collective rotation or a condensed vortex state, with no alternative. The 
ground state can undergo a kind of quantum phase transition when the shape 
of the interaction potential is smoothly varied. 
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Newly developed techniques [|l| for manipulating the strength of the effective interaction 
between trapped atoms 0] open the possibility to experimentally realize a Bose gas with 
weak interatomic interactions that are either attractive or repulsive. The limit {V) //icj^iO, 
where (y) is a typical expectation value of the interaction and hu is the quantum energy 
of the confining potential, is now reached experimentally via the Feshbach resonance Of 
special interest are the ground states of a rotating system P^jlj^l at given angular momenta 
L, the so-called yrast states ||5|,|^,f7i^,P,[T0| . 

So far, the problem was usually studied in dimension k=2 for the contact interaction 
y~a*^5(r), either attractive (a*^<0) or repulsive (a'*^>0) where a**'^ is the scattering length. 
The case a*'^<0 was solved analytically by Wilkin et al. [Q, while the case a*'^>0 was studied 
numerically by Bertsch and Papenbrock [§. The conjecture for the ground state wave 
function was confirmed analytically in |^, and generalized to a universality class of 
predominantly repulsive interactions A universality class of predominantly attractive 



interactions has been constructed in ||T0| 



Since experiments are now feasible [Q] in the weak-interaction limit close to the region 
where the scattering length changes sign, it becomes of importance to generalize these re- 
sults to arbitrary interactions. In this Letter the problem is solved exactly in k dimensions 
for arbitrary two-body central forces V{r). The final result is the following. In the func- 
tional space {V} of all possible interactions V{r), we may restrict our attention to those of 
physical interest, {V^hys}, of which we require that the force — ^ changes sign only once 
from repulsive at short to attractive at long distance: 

dV dV 

— < 0, r < R; > 0, r > R; R < 1. 1 

dr dr 

Since the crossover at R occurs for atomic reasons, we may assume it to be smaller than 
the trapping size, Rkuj^^^"^ (or R<1 in units uj=\). In addition, we assume that the force 



does not grow at r— ^oo and the quantum eigenvalue problem is free of divergencies [|rT|. It 
will be shown that the entire functional space {Vphys} is divided into two distinct classes 
{^phys\ ^"^^ {^phys\ '^^ (effectively) attractive or repulsive interactions. (The meaning of 
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'effective' in this context will be detailed below.) Within each class the energies of the 
yrast states depend in a simple way on the interaction while their wave functions remain 
the same. The ground states in the two classes differ qualitatively. When the interaction is 
changed from effectively attractive to repulsive, the transition between the two classes can be 
visualized as a quantum phase transition, with the relative interparticle angular momentum 
as a discrete order parameter. These exact analytical results are exemplified by analysis of 
Morse potentials with variable scattering length. 

The Hamiltonian of N spinless bosons of mass m in a /c-dimensional symmetric harmonic 
trap {k=2,3,...) reads 

H^m,.^).ivir.,^H,.V, (2) 

where Ho describes noninteracting particles, rj(pj) is the fc-dimensional position (momen- 
tum) vector of i-th boson, and V is the two-body interaction with rjj = |rj— rj|. Hereafter, 
\0l) denotes the ground state at fixed L and with maximum value of the conserved compo- 



nent, Lxy |12|. Other rotationally degenerate wave functions can be obtained from \0i) by 
applying the standard angular momentum ladder operators (e.g., L_=Lx—iLy, in k=3). In 
the limit hu^V, the determination of \0l) requires the diagonalization of the interaction 
within the space 

N 



Sz[^4\.J^\0), J2^^ = L, (3) 

i=0 

where Zi=Xi+iyi, S is the symmetrization operator, and \0)=e~^^^, in the convention 
h=m=uj=l. The states involving z*=Xi—iyi are separated from (^ by an energy 6E=nhu 
and their admixtures can be neglected for V<tihLJ. The dimensionality of the basis (0) grows 
exponentially with L 0. 

Within the subspace (||) the Hamiltonian (|^) is 

H = L + {Nk)/2 + W, (4) 

where the first (constant) terms come from Hq and W is the interaction V, projected 
onto the subspace (^. Using the integral transform 



w{l) 

we can write W in the form 



W = S Y.wiiij)S, kj = (at - a+)(ai - aj)l2 (6) 



i>j 



with lij the relative angular momentum between atoms i and j written in terms of the ladder 
operators o,t=^~^ ^^"^ ^i~'2~^'ir- "^^^ total internal angular momentum J=J2i>jhj is 
an exactly conserved quantity with eigenvalues J=^, j=0,2,3,..,L 0, commuting with the 
total angular momentum L=J2iO't(ii- 

We split [y] the projected interaction W into W = Wq + W s such that the first term 
is simple enough to find its lowest eigenvalue Eq and its associated eigenstate |0), 

Wo\0)=Eo\0). (7) 

The state |0) will also be the ground state of the total interaction iy=iyo+W^S' if (i) Ws is 
non-negative definite, 

Ws = Sj2vsitj)S>0, (8) 

i>j 

and (ii) |0) is annihilated by Ws 0, 

Ws\0) = 0. (9) 
In general, the operator Wq can be written as 

Wo = SY. Vo{h,)S, voil) = E c„r, (10) 

i>j m=0 

where Cm are hitherto unknown coefficients that need to be fixed such that conditions 
(0); (il)) (ID are satisfied. The condition (H) can be simphfied in the following way. 
The operator vs{ij) can be diagonalized in the unsymmetrized monomial states in Eq.(^ 
|m) = z'iZ2...z''^\0) via the substitution z^^ ^'^^ , ^j^ ^'^^ ■ Its eigenvalues are w{l)—Vo{l), 
Z=0,1,2,...,L ll^. Odd-/ eigenvectors are annihilated by S. Therefore, our strategy is to 



choose the Cm in such a way that some of the even-/ eigenvalues Xn=w {2n) —vo{2n) are zero 
and all others positive. This will guarantee that Ws is a. non-negative definite operator [Q, 
even though we cannot evaluate its eigenvalues directly. The candidates for the state |0) 
satisfying (|^) and (|^) can be sought among linear combinations of states S'|m)'s which obey 
f 5'|m)=0. 

The following solution for the Cm satisfies (||), (P): 

Co=w{0), C, = -A29{A2)/2, Cm>2=0, (11) 

where 9{x) = {0{x < 0), l(x > 0)} is the step function and A2n='w{0)—w{2n). The eigen- 
values (for A2<0) of the corresponding operator vs{ij) in (|]) read 

Xn = n9{A2)A2-A2n>0. (12) 

As will be shown below, the inequalities are satisfied for all interactions Vphys defined in (p. 
The ground state and its energy as a function of L {L<N) are 

L N 

\0,) = e-'^^^^SllT.i^^-^A^2)), (13) 

j=i i=i 

Nk NiN~l)wo e^NL A2e{A2) 
tjQ{L) = L + — - + 



2 2 4 

Here eL=6{L—l). The ground state depends on the interaction V{r) only via the sign of 
the control parameter 

2l~fc/2 ^oo 9/^4 



A, 



I nr)r'-V- *. (14) 



T(k/2) 

Among the potentials {Vphys}, the equation 

A2 = (15) 

selects the submanifold {V^*/ij,s} which we call separatrix. The se^^aratra divides the interac- 
tions Vphys{r) into the two classes {Vp^ys} and {V^y^}, differing by the form of the ground 
state: If A2<0, the ground state (|T3[) becomes 



|0-)=e-S^Z^ hj\Ol) = 0, (16) 



where Z=J2f§ [||] while for A2>0 and L>1 (||) gives 



|0+) = e-^^S-z,-z,...-z,, J|0+) = — 10+). (17) 

with Zj=Z—Zj. The internal angular momentum J can be viewed as a "discrete order 
parameter", zero in the domain of potentials V{r) with A2<0 and positive in the domain 



A2>0. For A27^0, (^) and (|l^) are non-degenerate ground states; within the separatrix, 
they are degenerate [ITsH . 

For k=2, we have rf=\z\f in Eqs.(|l^,|l^) and the |0j) accord with the results obtained 
for the attractive and repulsive [|], contact interactions, respectively. We arrive at an 
interesting conclusion: exact solution for the problem (|]) for arbitrary interaction is given 
by simple generalization of the results for the universality classes of predominantly attractive 
||T0| and repulsive interactions. 

By we denote the normalized occupancies z/„ {J2n^n=^) of the single-particle oscillator 
states z"- in the ground states |0^). The two behave qualitatively differently. For L=0 we 
have 1/^=1'^ =6n,o, that is the trivial (noninteracting) condensate. For 1<L< N we have 

L\[q{a - l)!2Fo(a,7;si) -FtQ!!2Fo(a + l,7;si)] 



sf Ar2"+%!(-7)!(Ar - L)\2Fo{-L, N-L + l;so)' 

Here 2-^0 is the hypergeometric sum [|l^, q=a+NP^Si and t=Si{2P+l) — l with 
Si={i+N)/ N'^ , a=N+'j, '-f=n—L, and (3=N —L—Nn. As L grows, the second excited z"^ 
and the first excited z^ states become populated (see Fig.l.). For Ly^O,N, the Uq , and 
V2 remain macroscopic [0(1)] in the limit A^^l. For L-^N the bosons condense in the 
state z^, forming a vortex. We call \0i) condensed vortex state In contrast, the state 
|0^) gives a binomial distribution z/^ = (^'!ln)\n\{2N-i)'- Peaked at n=^ p. We call |0^) a 
collective rotation state: its angular momentum comes from the collective factor Z^. An 
interesting way to view such states as rotated noninteracting condensates was presented in 



The control parameter A2 measures the balance between repulsion and attraction in 
V{r). For zero-range interactions, 2^-1% ^2 coincides with the scattering length in Born 
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approximation, ^1^^= J^V{r)r^~^dr. For finite-range interactions it differs from af,'^; and it 
can be considered as a modified scattering length whicli separates tlie two universality classes 
of effectively attractive and effectively repulsive interactions. 

For A2<0, the interaction contribution to the ground states energy (^) is negative and 
L- independent. It grows with A2. At the critical point A2=0, the derivative dEQ{L)/dA2 
has a jump '^^"^ and Eq(L) becomes L-dependent for A2>0. In Fig.2, we show the results 
for the Morse potential VMir)= e ' ' —2e^^. We display Vm for two sets of parameters 
a and R, and the resulting /-even part, We, of the transform (|[) 

w{l) = r(/ + — [e - dl^ilk - 2e rf^^;4J , (18) 

where di^^^=e^^^^Di{s) with -D/(s) the parabolic cylinder function ||16[. It is seen that 



Xn=w{2n)—VQ{2n)>0 in both cases. When attraction prevails, A2<0, We has a minimum 
at /=0 and the system resides in its lowest-energy state |0^) for a minimum possible value 
of J=Y^hj{=0). When repulsion starts to prevail, the quantity A2 becomes positive, the 
minimum of We is shifted to Z>0, and the system finds a lower energy for a non-zero value 
of J in the ground state |0^). The behavior of We{l) at the critical point A2=0 resembles 
that of a thermodynamic potential in a second-order phase transition |T^. In contrast to 



thermodynamical systems, this "phase transition" is insensitive to k and it happens for finite 
A^. Panel c) shows the bifurcation point in Eq{L). 

Fig.2 d) shows the phase diagram for the Morse potentials in the parametric space {a,R). 
The separatrix (|T5|) between {V^} and {V^} reduces to a curve defined by the following 
relation between a and R 

Ria) = a ln(0itV0^t\), ^ - ^41- (19) 

Using Eqs. (p!8|JT9| ) , Ws can be shown to be non-negative definite (0) throughout: For values 
of R<R{a), we obtain A2<0 and A^>0; for R>R{a), we have A2>0 and A^>0. Thus all 
Morse potentials, whatever their values of a and R are covered by (|T3]) and fall into two 
distinct classes depending on whether A2 is positive or negative. [Note that for R>1, the 



Vm lies outside the class (|T|)!] We obtained similar results for other two-parametric potential 
families. 

Is this situation generic? The full set of physical interactions (|1]) forms a dense functional 
manifold. In general, it cannot be described by a finite number of parameters. Let {V} 
be the complete manifold of all potentials [not necessarily the physical potentials (|I])] as 
indicated by the big disc in Fig. 3. Within this extended manifold, we can still define the 
subclasses that have A~>0 and A^>0 with the ground states (|16D and ([T7|), respectively. 
Their boundaries A~ and can in general be distinct, leaving room marked by when 
the ground state is not (p^Sf). We will show that there is no such gap for all the physical 
interactions (|l|), confirming that the result (|13]) covers the whole universality class {Vphys} 



and that the separatrix A2=0 ([T5| ) divides the interactions {Vphys} (small disc in Fig.3.) into 
the two classes, with the ground states ([T6|) and (0), respectively, with no other alternative. 
The proof requires to show that for any potential in {Vphys} A~>0 for A2<0 and A^>0 for 
A2>0, Eq.(Q). 

The quantity A2„= J^(f)2nFdt is the work done by the modified force F(j)2n=-^^^^^(p2n 
to separate a pair of particles. Here 02n=Z]m=i ''r(fc/2+m) positive definite. By (|T]) the 
work is the sum of the positive and negative areas 



rT roc ^2 

A2n = / hn(j)2Fdt + / hn<p2Fdt, r = —. (20) 
Jo Jt 2 

where hn=^ are positive-valued functions. If A2<0, the second (negative) term prevails for 
n=l. It will then prevail for any n>l, because the functions /i„>i increase monotonically. 
That is, if A2<0, then all A2„<0 and A->0. If A2>0, the first (positive) term in (|0D 
prevails for n=l. We introduce the monotonically decreasing functions h=n—hn- They are 
positive for t<t„ and negative for t>t„ with all tn>l. By ([1|) we have T<tn and we write 

A+ = / K(t)2Fdt + / K(t)2Fdt + / K(t)2Fdt. 

Jo Jt Jtn 

The second term is the only negative contribution. The corresponding area is however 
smaller than that of the first term, because hn decreases monotonically and T<tn- Thus 
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A^>0 for A2>0. To sum up, the inequalities (0) hold throughout, and the solution ([I3D 
covers the universality class {Vphys}- 

By similar arguments, one can append the class {Vphys{r)} ([1|) by potentials with con- 
stant sign of F and '^'^if'\ like 6{r), i, log{r), e"*^/", e"'''/"' etc. 

In summary, we have solved the problem of rotating ground states of weakly interacting 
trapped Bose atoms. The resulting phase diagram gives a complete classification of these 
ground states. The conclusion is that the ground state is either a collective rotation or 
a condensed vortex state; in both cases analytical expressions are known for the ground- 
state wave function and its energy and these are valid whatever the form of interatomic 
interaction, as long as it has reasonable properties. These results can be of interest in view 
of recent experiments with trapped atoms and can be applied to other systems. 

The work was supported by CEA (France) and FAPESP (Brazil). 
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Figure Captions 

FIG. 1. Left: Asymptotic values of as functions of L/N for r;,=0,l,2 at A^^l. Right: 
Vn as functions of n for L=12 (dash) and for L=24 (sohd) for A2>0 and A2<0 (A^=30). 

FIG. 2. Morse potential V Mij) (a) and Weil) from ([T8|) (b) for the two sets of parameters: 
"l"(a=0.8,i?=1.0,A2>0) and "2" (a=1.15,/2=1.55,A2>0), A;=3. Symbols connected by solid 
lines show We- The Vq (|I^,|n|) are shown by dashed lines, (c) The interaction contribution 
to the ground state energy (|T^, as function of R (for a=l/3, A;=3 and A^=4). (d) Phase 
diagram for family of the Morse potentials on the a— i?-plane. The separatrices (|TUp for k=Z 
(solid) and k=2 (dash). 

FIG. 3. Schematic phase diagram in functional space {V{r)}. Within {Vphys}, the 
boundaries A~ and A"*" merge with the separatrix A2=0 ([T5|) . 
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FIG. 1. FIG. 1. Left: Asymptotic values of z/+ as functions of L/N for ra=0,l,2 at A^^l. 
Right: Vn as functions of n for L=12 (dash) and for L=24 (solid) for A2>0 and A2<0 (A^=30). 
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FIG. 2. Morse potential VmI?") (a) and We{l) from (18) (b) for the two sets of parameters: 
"l"(a=0.8,i?=1.0,A2<0) and "2"(a=1.15,i?=1.55,A2>0), k=3. Symbols connected by solid 
lines show We- The Vq (10,11) are shown by dashed lines, (c) The interaction contribution 
to the ground state energy (17), as function of R (for a=l/3, k=3 and A^=4). (d) Phase 
diagram for family of the Morse potentials on the a— i?-plane. The separatrices (19) for k=3 
(solid) and k=2 (dash). 
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FIG. 3. Schematic phase diagram in functional space Within {Vphys}^ the 

boundaries A~ and A"*" merge with the separatrix A2=0 (15). 
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